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ABSTRACT 

Using a high cadence imaging spectropolarimetric observation of a sunspot 
and its surroundings in magnetically sensitive (Fe I 6173 A ) and insensitive (Fe 
I 7090 A ) upper photospheric absorption lines, we map the instantaneous wave 
phases and helioseismic travel times as a function of observation height and in- 
clination of magnetic field to the vertical. We confirm the magnetic inclination 
angle dependent transmission of incident acoustic waves into upward propagat- 
ing waves, and derive (1) proof that helioseismic travel times receive direction 
dependent contributions from such waves and hence cause errors in conventional 
flow inferences, (2) evidences for acoustic wave sources beneath the umbral pho- 
tosphere, and (3) significant differences in travel times measured from the chosen 
magnetically sensitive and insensitive spectral lines. 

Subject headings: radiative transfer — Sun: helioseismology — Sun: oscillations 
— Sun: surface magnetism — sunspots 



1. Introduction 



Accounting for directly observable photospheric wave evolution within a sunspot while 
inferring subsurf ace conditions is an important and challenging problem in the seismol- 
ogy of sunspots (IBogdan et al.lll998l ). The characteristic association between non- vertical 
magnetic fields and acoustic wave propagation at frequencies well below the photospheric 
acoustic cutoff of ~ 5.2 mHz, observed in sunspot penumbrae (IMcIntosh fe Jefferiesl 12006 



Rajaguru et al.l 120071) as well as in disparate wave dynamical phenomena in the solar at 



mosphere foe Pontieu et al.l 12004 120051 : Ijefferies et al.l 120061 : Ide Wiin et al.l l2009h ha s fur- 



ther highlighted this problem. Several indirect influences, due to physical (jWoodardl 119971 ; 
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Bogdan et al. Il998h as well as analysis specific re asons ( jRajaguru et al.|[2006l IParchevsky &: Kosovichevl 
20071 ). of p-mode absorption (IBraun et al.lll987l ) per se have been shown to manifest as ap- 
parent flow and wave speed signals in local helioseismic measurements; however, there have 
been none studying the direct contribut ions of propagating waves. The "inclin ed magnetic 
field effect" in helioseismic signatures ( jSchunker et al.l 120051 ; IZhao et al.l 120061 ) has indeed 
been shown to arise from magnetic-field-aligned wave motion caused by an incident acoustic 
wave, but the interpretations relied on a viewing-angle-dependent geometric relation between 
the magnetic field and wave motion, whose phase shifts are independent of those arising from 
wave progression in height. Viewing-angle-dependent changes in observation height too could 
cause different wave phases at different positions w ithin a sunspot and po tentially could be 
misinterpreted as the above effect, as suggested by lRajaguru et al.l ( 120061 ). 



In this Letter, we derive explicit observational proofs for helioseismic contributions from 
propagating waves within a sunspot and from acoustic sources located beneath its umbral 
photosphere. We also show significant differences in travel times measured using velocity 
data from magnetic and non-magnetic lines. 



2. Observations and Analysis Methods 



We performed imaging spectropolarimetry using the Interferometric Bl-dimensional 
Spectrometer (IBIS) installed at the Dunn Solar Telescope of the National Solar Observatory, 
Sac Peak, New Mexico, USA. IBIS has spectral and spatial resolutions of 25 mA and 0".165, 
respectively, and has a 80" diameter (^60 Mm) circular field of view (FOV). We observed 
a medium sized sunspot (NOAA AR10960, diameter rs 18 Mm) located close to the disk 
center (S07W17) on 2007 June 8. Our observations involved scanning and imaging in all 
the Stokes profiles (J, Q, U, V) of magnetic Fe I 6173.34 A and in Stokes / of non-magnetic 
Fe I 7090.4 A , with a cadence of 47.5 s. A 7 hr continuous observation was chosen for our 
analysis. The spectral images were dark subtracted, flat-fielded, and re-registered to remove 
atmospheric distortions, which were derived from white-light images recorded simultane- 
ously. The spectral calibration consisted of removing a quadratic center-to-edge wavelength 
dependence, and correcting for the transmission profiles of the prefilters. Polarization cal- 
ibration of the magnetic line data was based on the telescope and instrument polarization 
matrices (jCavaUinil 120061 ) . 



Similar to iRajaguru et al.l (120071 ). we extract line-of-sight (LOS) velocities of plasma 
motions within the line forming layers from the Doppler shifts of line bisectors. We use 10 
bisector levels with equal spacing in line intensity, ordered from the line core (level 0) to the 
wings (level 9), and derive 10 velocity data cubes, Vi(x, y, t)(i = 0, 9), for each line. For the 
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magnetic line, we use the average o f bisector velocities from the left (I + V) and right (I — V) 



circu lar polarization (CP) profiles (ISankarasubramanian fe Rimmeld 120021 ; Idel Toro Iniesta 



20031 ) and those from the I profile for the non-magnetic line. The 10 bisector levels span 



the height range within the line formation region in an unique one-to-one way. Based on 
the Maltby-M umbral model atmosphere ( Malt by et al.l Il986 ). the Fe I 6173.34 A line 
formation is reported to span a heig ht range of 20 km ( wings) to 270 km (line core) above 
continuum optical depth r r = l level ([Norton et al.l 120061 ) . and a very similar range for Fe I 



7090.4 A (jStraus et al.l 120081 ). We choose about 1 hr long observation from the best seeing 
interval (the first 3 hr), and do Milne-Eddington (M-E) inversions of the (temporal) average 
of Stokes profiles of the mag netic line to obtain magn etic field B(x,y), its LOS inclination 
j(x, y) and azimuth ip(x,y) (jSkumanich fe Liteslll987l ). Because of low polarization signals, 
the above inverted quantities are noisy outside of the sunspot (see Figure 1), and we cut out 
only the spot region for use in our analysis. 



2.1. Instantaneous Wave Phases and Helioseismic Travel Times 



Instantaneous wave phases in the form of phase shifts 50j i o(^)=Phase[Vj(z/)VQ(^)], 
where v is the cyclic frequency of a wave and V is the Fourier transform of v, due to wave 
progression between two heights corresponding to any one of the bisecto r levels i = 1, 2, 9 
and level (the top most layer) are calculated (IRajaguru et al.l 120071 ) . Since we want to 
study mainly the p-modes and compare 5<p with helioseismic travel times, we filter out the 
/-mode. We take median values of 5<p over the p-mode band (2-5 mHz) or over bands 
of 1 mHz full width at half-maximum (FWHM) centered at every 0.25 mHz (to study any 
frequency dependence). Signals over space are studied using 7 or LOS magnetic field B L qs- 
For this work we focus on studying the 7 dependence of 50 and use 3° bins in 7 (see Section 
O below). 

The 10 different data cubes from each line are run through a standard p-mode time- 
distance analysis procedure in center-annulus geometry (IRajaguru et al.ll2004f ) . In addition, 
we apply frequency filters the same way as for 0. Trav el time maps are calcula ted for five 
travel distances A = 6.2, 8.7, 11.6, 16.95, and 24 Mm (see lCouvidat fe Birchl J2009h for details 
on phase speed filters). For this Letter, we focus on analyzing the results for A = 16.95 Mm, 
because, given the sizes of observed region (radius ~ 29 Mm) and the spot (radius ~ 9 Mm), 
this is the optimum A that facilitates distinguishing clearly the ingoing and outgoing waves in 
the sense of their interactions with the spot. We also perform a double-skip annulus-annulus 
geometry measurement, for A = 16.95 Mm, which avoids use of oscillation signals within the 
sunspot, for the diagnostic checks presented in Section 4. Height dependent contributions 
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to outgoing and ingoing phase travel times r + and t~ from within the line forming layers 
are determined using St^ — — (i — 1, 9). To facilitate comparisons with 5<p, we 
average St^ too over 3° bins of 7. All pixels outside of 10 Mm radius centered on the spot, 
for convenience (in Figures 2-4), are assigned a 7 value of 90°; all wave quantities averaged 
over these set of pixels are taken as that of 'quiet-Sun', and a subscript q is used, where 
necessary, to identify them explicitly. 

2.2. Error Analysis 

The oscillation signals are inherently stochastic due to such nature of acoustic sources 
and that of the background medium. A typical measurement of a wave quantity, hence, here 
either in S<fi or St, carries a random error. Assuming that all pixels with the same 7 or within 
a small range of 7 comprise independent measurements of the same S<p and St, we take the 
mean over these pixels as our best estimate and study its variation against 7. Error estimates 
for 7 from M-E inversions fall in the range of 0°.85 - 2°, with mean values of 1.5° over the 
umbra and 1° over the penumbra. A bin size of 3° in 7, referred to in previous subsection, 
is found to be optimal to accumulate a statistical sample of measurements while being small 
enough to not bias them through their variation against 7 itself. So, error estimates for the 
means S(f> or St are their standard errors given by cr/yjn, where a is the standard deviation 
in n number of measurements (i.e. pixels falling within a given bin in 7). 

3. Origin and Seismic Contributions of Propagating Waves 

Oscillation signals in photospheric Doppler velocities, in general, would consist of evanes- 
cent waves caused by the p-modes trapped below it and various propagating waves traveling 
in different directions. The <50i,o arise only due to the latter propagating waves. The surface- 
s-mode) and atmospheric-gravity waves also show vertical phase propagation in the pho- 
tospheric layers, but we have filtered them out here in our analysis. In general, Scf)^ should 
receive contributions from waves locally generated and those generated elsewhere (e.g., from 
the quiet-Sun) but get "converted" by the magnetic field to propagate upwards upon inci- 
dent on it. On the other hand, the ingoing travel times St~ get contributions solely from 
the latter helioseismic ones, which in our current analysis case are those traveling from a 
distance of A = 16.95 Mm in regions surrounding the spot. We show in Figure 2 S<ps,o and 
St$ , due to wave evolution within the region bounded by the wing (level 8) and core (level 
0) formation heights, against 7. The v values marked in the panels of Figure 2 are the central 
frequencies of 1 mHz band filters used. Keeping in mind that <50s,o have contributions from 
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a larger set of waves (as discussed above), results in Figure 2(a) for the magnetic line show 
a surprising amount of correlation between the two measurements, and moreover exhibit a 
strikingly similar 7 dependence. These results immediately reveal several interesting aspects 
of magnetic field - acoustic wave interactions: (1) first of all they confirm that helioseismic 
waves incident on the sunspot see themselves through to higher layers of its atmosphere 
with a striking dependence on 7: a coherent transmission of incident waves happen, peaking 
around 7 ~ 30°, maintaining a smooth evolution of time-distance correlations; (2) remem- 
bering that CP profiles of the magnetic line have maximum sensitivities for velocities within 
vertical magnetic field, it is seen that a large fraction of waves propagating upward within 
such field are due to helioseismic waves originating at distant locations; and, (3) provide di- 
rect evidences that ingoing wave travel times would cause observing height dependent signals 
in flow inferences from travel time differences. The non- magnetic line (Figure 2(b)) yields 
very little correlations between <50g,o an d St$ ; however, the helioseismic measurements <5t 8 _ 
agree well with that in Figure2(a), except at high v, thus reinforcing inferences (1) and (3). 
We speculate that there are substantial wave motions, locally generated, perhaps within the 
non-magnetic gaps or weakly magnetized penumbral region, whose signatures are missed 
in the CP profiles of the magnetic line. At v > 4 mHz, helioseismic signatures within the 
sunspot get markedly different in magnetic and non-magnetic lines. 

To affirm the reader that the signals analyzed in Figure 2 (as well as other Figures) 
are due to the height evolution of wave phases and not due to any other wave correlations, 
random or spurious, in Figure 5 we show maps of 5(f) i;0 and 5t~ for the full p-mode band 
(2-5 mHz), whose variations against i = 1, 2, 8 are in the x-direction of the maps. The 
gradual increase in the values of £0i,o and 8t~ with that in height separation is obvious in 
this figure. We further note that halving the bin size in 7 to 1.5° changes negligibly the 
mean values 5<p and 5t and their variation against 7 studied in this work, but increases their 
error estimates by roughly times. 



4. Wave Sources Beneath Umbral Photosphere 

Outgoing waves at a given measurement location, in general, would consist of those 
locally generated and those generated elsewhere undergoing reflection at the photosphere 
directly below it. These latter component would be seen in neither 50^0 nor ^toi as they are 
evanescent at the observing height. For locally generated waves, circular wavefronts from 
a source, while their upward propagating parts see themselves up through the magnetic 
field, would cause outgoing wave correlations yielding distinct signatures in 5t^ (see Figure 
3(b)). Results in Figure 3(a), for 5t^ from both the magnetic and non-magnetic lines, do 



-6 - 



indeed provide such a diagnostic: outgoing waves starting at higher height (line core) within 
the sunspot atmosphere and reaching the quiet-Sun at the chosen A have shorter travel 
times than those starting at a lower height (line wings) and reaching the same quiet-Sun 
location; since this is simply not possible, the only explanation for this observation is the one 
contained in our previous sentence and illustrated in Figure 3(b), viz., outgoing wave time- 
distance correlations are predominantly due to waves directly from sources just beneath the 
sunspot photosphere when oscillations observed within it are used. We also note the different 
7 dependence of 5r + , as compared to that of 5r~ , possibly due to circular wavefronts from 
local sources having a range of incident angles with the magnetic field. 

To confirm the above and to check the extent of contributions from magnetic patches 
surrounding the spot, we measured travel times in double-skip annulus-annulus geometry 
that avoids waves observed within the sunspot. Here, diametrically opposite points on 
the annulus are correlated and azimuthally averaged; hence only the mean travel times are 
measured and assigned to the center points, directly beneath which the waves reflect between 
their two skips. The results in Figure 4 show that St$ s are indeed small, compared to St$ . 
Hence, the small scale magnetic patches surrounding the spot do not contribute much to St^q 
and hence to St^ too, validating our inferences above and in Section 3. As to magnitudes 
of St (Figures 3a and 4), it is interesting to note that they are about thrice those of St^ 
(Figure 2) (in the main p-mode band of 2 - 4 mHz). This difference could arise from non- 
circular expansion of wavefronts possibly due to two causes: (1) the differences in physical 
conditions, due to material flows as well as sound speed, in the wave path regions H c -H„, 
and H^,-H^ (refer to Figure 3(b)), and (2) the motion of the sources beneath the umbral 
photosphere. In the lower panel of Figure 4, we compare mean travel time perturbations 
measured from line core and wing bisector velocities, t™-t!^ and t^-tI^, with those of half 
the double-skip travel times, {tI s -t%1)/2, from wing bisector velocities. Interestingly, almost 
all the differences between double- and single-skip travel times appear to come from the line 
formation layers, and hence are observation height dependent. These results also show that 
height difference between the source location and the wave reflection layer beneath the spot 
are very small. 



5. Discussions and conclusion 

Almost all time-distance helioseismic analyses proceed under the working assumption 
that wave signals at observation heights are evanescent and hence oppositely directed wave 
paths involving photospheric reflections at two separated points are of identical path length. 
This assumption is basic to the inferences on flows and wave speed from travel time differences 
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and mean, respectively. In an early theor e tical study, accompa nied by attempts to model the 
helioseismic observations of lBraunl ( 119971 ). iBogdan et al.l ( 119981 ) showed the influences of both 
the p-mode forcing of, and spontaneous emissions by, sunspots on acoustic wave travel times. 
Our analyses here have yielded transparent observational proofs for both effects, for the first 
time, with important new perspectives: (1) the process of transformation of incident acoustic 
waves into propagating (magneto)-acoustic waves up through the magnetic field happen 
in a coherent manner allowing a smooth evolution of time-d istance correlations and, in 
agreement with several rece nt theoretical and numerical studies ( lCalryll2005l ; ICrouch fc Cally 
20051 ; ISchunker et al.l 120061 ). this process depends on the inclination angle (7) of magnetic 
field to the vertical, and (2) outgoing waves from acoustic sources located just beneath the 
sunspot photosphere add important additional contributions for both mean travel times and 
differences. Our results have also shown observational prospects for consistently accounting 
for the above effects in sunspot seismology, viz. the indispensability of imaging spectroscopy 
to extract wave fields so as to be able to correctly account for the wave evolution within 
the directly observable layers of sunspot atmosphere. Current limitations in making such 
observations over large enough FOV do not allow us to perform seismic inversions reliably. 
However, the analysis methods followed here point ways to a consistent and much improved 
observational determinations of structure and flows beneath sunspots once our instrumental 
capabilities improve. These observational avenues also promise a close scrutiny of various 
theoretical ideas and models of acoustic wave - magnetic field interactions and those of the 
associated MHD waves and their propagation characteristics. 
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Fig. 1. — Maps of continuum intensity I c , and magnetic field B, its azimuth ip, and LOS 
inclination 7 over the sunspot region from M-E inversions of Fe I 6173 A Stokes profiles. 
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a) Fe I 6173 A, CP Wings vs Core b) Fe I 7090 A, I Wings vs Core 




Fig. 2. — Instantaneous phase shifts, 5(psfl{^) (in red), and changes in ingoing wave travel 
times, 5t$ (in green), due to wave propagation between the formation heights of wings (20 
km) and core (270 km) of Fe I 6173 A {panel a), and of Fe I 7090 A {panel b) against 7 of B. 
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Fig. 3. — {Panel a) Changes in outgoing wave travel times, 5r^ , due to wave propagation 
between the formation heights of the wings (20 km) and the core (270 km) of Fe I 6173 A (in 
red) and 7090 A (in green) lines as a function of 7 of B. {Panel b): an illustration depicting 
wavefronts from acoustic sources beneath the umbra, wave paths and line formation heights 
(see text for details). 
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Fig. 4. — Upper panel: comparing outgoing wave travel times 5t$ and double-skip mean 
travel times 5t$ s from within the line formation region. Lower panel: changes, with respect 
to quiet-Sun, in half of double-skip travel times (violet) and in mean travel times at formation 
heights of line wings (green) and core (red). 
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Fig. 5. — Maps of <50 iiO (top row) and 5t~ (bottom row) for the full p-mode band (2-5 mHz). 
The variation in the x-direction of the above maps is for the bisector levels % = 1,2,..., 8 
corresponding to progressively deeper locations in the sunspot atmosphere with respect to 
% — 0, which is for the line core forming at 270 km above the photosphere. 



